Dimer motion on a periodic substrate: Spontaneous symmetry breaking and absolute 

negative mobility 
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We consider two coupled particles moving along a periodic substrate potential with negligible 
inertia effects (overdamped limit). Even when the particles are identical and the substrate spatially 
symmetric, a sinusoidal external driving of appropriate amplitude and frequency may lead to spon- 
taneous symmetry breaking in the form of a permanent directed motion of the dimer. Thermal 
noise restores ergodicity and thus zero net velocity, but entails arbitrarily fast diffusion of the dimer 
for sufficiently weak noise. Moreover, upon application of a static bias force, the dimer exhibits a 
motion opposite to that force (absolute negative mobility). The key requirement for all these effects 
is a non-convex interaction potential of the two particles. 

PACS numbers: 05.45.-a,05.40.-a,05.60.-k 
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I. INTRODUCTION 

The dynamics of a dimer in the presence of a periodic 
potential is of interest in a variety of different contexts. A 
first example is a diatomic molecule adsorbed and mov- 
ing on a crystal surface JjKa] , e.g. during crystal growth 
or thin film formation 0t3- In many cases, the dimer 
is able to explore the entire two-dimensional surface, but 
in some cases also an effectively one-dimensional motion 
arises, e.g. due to a strong crystalinc anisotropy, apro- 
nounced step, or a channeled (110) metal surface [a, 0j- 
As an example, Fig. 1 illustrates a Si2 or Gc2 dimer con- 
fined between the rows of the reconstructed Si(001) or 
Ge(001) surface, see [1] and references therein. Motion 
of two coupled atoms (or defects, nanoparticles, etc.) in 
a three-dimensional (atomic, optical etc.) lattice is yet 
another option along these lines. Further important ex- 
amples are simple models of processive molecular mo- 
tors consisting of two elastically coupled "heads" (motor 
domains) and walking along a polymer filament [9l-ll8j. 
Artificially manufactured set-ups include colloidal dou- 
blets [lj| in periodic arrays of magnetic [2(| or optical 
[2lj traps. Last but not least, dimer models are well es- 
tablished to describe the sliding friction of two coupled 
asperities at the nanoscale |22h25| . but also arise in less 
"obvious" systems like superionic conductors [26[, disso- 
ciated dislocations 0, and SQUIDs [11, 111- 




FIG. 1: (Color online): Schematic illustration of the confined 
motion of a Si2 or Ge2 dimer (big red spheres) within the rows 
of a 2x1 reconstructed Si(001) or Ge(001) surface (small grey 
spheres) . 



In most of the previous theoretical studies, the dimer 
was driven out of equilibrium by some external forcing 
and the spatial inversion symmetry of the system was 
broken a priori in one or another way. Most promi- 
nently, in the context of so called ratchet effects, a spa- 
tial anisotropy of either the substrate potential or of the 
dimer itself was assumed @, El [lJ-Gj, H3-53 ■ Another 
important direction concerns the behavior of drift and 
diffusion when the spatial symmetry is broken by an ex- 
ternally applied, constant "bias" -force @, [H HJ5I 51 . 
Studies of fully symmetric dimer models have focused on 
the diffusive properties, since a systematic directed dimer 
motion seems - at first glance "obviously" - not possible 
in such systems 0, H, 0, 0, & HH SI ■ Exceptions are some 
idealized conservative (Hamiltonian) models 0, SH and 
the recent paper (52|, to which we come back in Sec. X. 

In our present study, the main focus is on dimers 
consisting of identical particles and on spatially peri- 
odic systems without any intrinsically broken symme- 
try. To avoid unnecessary complications, we furthermore 
assume that inertia effects are negligible (overdamped 

limit) n, iii umi [ms m m, isa m siii, and 

that the dimer is driven out of equilibrium by a sinu- 
soidally oscillating force (HHH, IS SI IM 5l SI SI . 

Our first main result will be that - in spite of the spatial 
symmetry of the system - the deterministic dynamics still 
admits a systematic net motion of the dimer via spon- 
taneous symmetry breaking. In particular, ergodicity is 
broken, and thus the direction of the spontaneous trans- 
port depends on the initial conditions. Upon including 
thermal noise, ergodicity is restored, and thus the aver- 
age dimer velocity must be zero. But now - quasi as a 
precursor of the spontaneous deterministic transport - 
the diffusion coefficient diverges for asymptotically weak 
noise. Finally, we will demonstrate that upon applica- 
tion of a static bias force, the system may respond with 
motion opposite to that force, i.e. it exhibits so-called 
absolute negative mobility. The main and indispensable 
prerequisite for all these effects is a non-convex interac- 
tion potential of the two particles. This quite unexpected 
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and remarkable condition also seems to explain why such 
effects have not been "accidentally" discovered in previ- 
ous studies of overdamped dimcr motion. An exception 
is the recent absolute negative mobility effect in Ref. (29| 
which will be discussed in more detail in Sec. X. 



II. ONE-DIMENSIONAL MODEL 

With the exception of Sec. IX, we will always focus 
on the one-dimensional motion of two coupled, identical 
particles with coordinates X{ (i = 1,2), modeled by the 
Langevin equations 

• /a dU t ot(xi{t),x 2 (t)) 
T]Xi{t) = — \-F + f{t)+€i{t). (1) 

In doing so, inertia terms m,Xi(t) are considered as negli- 
gible (overdamped limit), rj is the viscous friction (damp- 
ing) coefficient of the single particles, and F a static 
"bias"-force. Furthermore, f(t) represents an external 
driving with time period r, satisfying the symmetry con- 
dition 

/(f + r/2) = -/(*), (2) 

and thermal fluctuations are modeled by independent 
Gaussian noises £i(i) of zero mean and respecting the 
fluctuation-dissipation relation 

(&(*)&(*)) =2 V k B T5 ij S(t-a) , (3) 

where fee is Boltzmann's constant and T the ambient 
temperature. The potential energy Utot(%i, x 2 ) is of the 
form 

U t ot(xi,x 2 ) = U( Xl ) + U(x 2 ) + W(xi ~ x 2 ) (4) 

with a spatially symmetric and L-periodic single-particle 
potential U{x) and an interaction potential W{y) which 
will be further specified later on. 

For the sake of convenience, we will often adopt units 
of time, length, temperature, and energy so that 

rj = 1, L = 2ir, ks = 1, max U (x) — min U(x) = 2 . (5) 

X X 

When specific examples will be needed, we will further- 
more focus on 

U{x) = cos(x) (6) 
f(t) = A S m(ut) (7) 

with lo := 2ir/r, and either a Lennard- Jones or a quartic 
interaction potential of the form (see also Fig. 2) 

W(y) = (\/y) 12 -2(\/yf , (8) 
W(y) = c(j/-A) 3 (y-A-2) , (9) 

with positive parameters c, A, A. We expect and explic- 
itly verified in a few cases that a large variety of similar 
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FIG. 2: (Color online) (a) The Lennard- Jones interaction 
potential ((HJ for A = 1.25 L. (b) The quartic interaction po- 
tential © for c = 0.428 and A = 0.26 L (see also Fig. 0. 



potentials U(x), W(y) and forces f(t) lead to qualita- 
tively unchanged results (see also Sees. VI and IX). 

The Lennard-Jones potential is a well established 
prototype model for the interaction of atoms. In our 
present context, it has the only, and mainly puristic, 
shortcoming that for any finite thermal noise, a dimer ul- 
timately dissociates after sufficiently long time into two 
separate monomers. Our second example (J9j) is free of 
this problem and, moreover, exhibits a slightly richer 
"infra-well" structure of the potential landscape, see Fig. 
2. In particular, it has a metastable inflection point at 
y = A, thus mediating between a "purely" monostablc 
and a "genuine" bistable interaction potential. We thus 
can expect that the results we will find for © are still 
qualitatively recovered for suitable monostablc as well as 
bistable interaction potentials (see Sec. VI. C). 

Note that interaction potentials in real systems are ex- 
pected to satisfy the symmetry W(— y) = W(y). In con- 
trast, we also admit non-symmetric W(y), e.g. in (J9j> , for 
the mere sake of convenience later on. However, when 
dealing with asymmetric W(y), we tacitly restrict our- 
selves to cases for which a change of sign of x\ — x 2 in 
is a negligibly rare event in the dynamics ([1]), e.g. 
by admitting only sufficiently large c- values in In 
other words, the two monomers do not pass each other 
in practice, and thus only y- values of one sign, say y > 0, 
actually matter in W(y). 

Later on, our main focus will be on non-convex inter- 
action potentials W(y). Various physical realizations of 
such potentials will be discussed in Sec. V. Depending 
on the particular system, the external driving forces F 
and f(t) may be realized by hydrodynamic flows, trav- 
eling optical potentials, electrical fields, magnetic fields, 
or mechanical agitation. 



III. BASIC NOTIONS 

Given the model ©-(HI, ^ ^ s na tural to go over to 
center of mass and relative coordinates 

x := (x\ + x 2 )/2 , y := xi - x 2 . (10) 
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In the special case (0-([7]), their equations of motion take 
the form 

x(t) = sin(x)co8(y/2) + Asin(u)t) + F + Vr£x(t) (11) 
y(t) = 2cos(x)sin(y/2)-2W'(y) + 2Vfi 2 (t) (12) 

where £1,2 (0 arc independent, delta-correlated Gaussian 
noises. 

A first main quantity of interest is the time averaged 
velocity of the dimer's center of mass, 

t 1 r* 

v := — lim — / dt' x(t') , (13) 

L t~>00 t Jq 

expressed as a dimensionless multiple of the spatial and 
temporal periods L and r — 2it/uj for reasons of conve- 
nience later on. 

For any finite noise strength T, the dynamics ([I]) is 
ergodic, i.e. the velocity (TIB")) is (with probability one) 
equal to the ensemble averaged velocity. In particular, 
the velocity is independent of the initial conditions Xi (0) 
and of the realization of the noise £j (t) . 

In contrast, in the deterministic limit (T = 0), the 
velocity (fT3")) may - as we will see later - still depend on 
the initial conditions Xj(0), implying that ergodicity of 
the dynamics ([T]) is broken. 

A second quantity of interest is the diffusion coefficient 
of the dimer's center of mass, 

D := lim <* 2 (*)>-<°(*)> 2 , ( 14) 

t^oo 2t 

where (■) indicates an ensemble average over different re- 
alizations of the noise &(£). For any T > 0, ergodicity 
implies that an additional averaging over different ini- 
tial condition Xi(0) is possible but actually superfluous 
in 0. 

Next, we turn to the notion of spontaneous symmetry 
breaking. Generally speaking, this notion refers to a de- 
scription (model) of a system in terms of some mathemat- 
ical equations which exhibit a certain symmetry (invari- 
ance), while their solutions (system states) do not exhibit 
this symmetry. The simplest example is the equations 

a + b= 1 , ab = (15) 

which are clearly symmetric in a and 6, while their so- 
lutions, namely (a, b) = (1,0) and (a, 6) = (0,1), both 
break this symmetry. It also exemplifies the general rule 
that different symmetry breaking solutions arc related to 
each other via the original system symmetry. The sim- 
plest physical example is tossing a coin, illustrating that 
the decision about which state the system actually se- 
lects is usually rooted in the initial conditions and/or 
some "random effects" . 

Returning to our dynamics (H])-((4|), we observe that it 
satisfies the following basic symmetry property: If £1,2 (£) 
is a solution of (JTJ) then - due to the symmetry of f(t) 
in © and of U(x) mentioned below ((U) - also x\ 2 (i) : = 



—3>2,i(t + t/2) will be a solution of (Q}, provided F is 
replaced by -F and f 1)a (t) by f} 3 (t) := -£2,1$ + t/2). 

For T — (deterministic limit) and F = (unbiased 
system) the dynamics ©-(HI) is thus - up to an irrelevant 
shift of the time-origin and a relabeling of the (identical) 
particles - symmetric under spatial inversion. It follows 
that for any possible value of the velocity v according to 
([TBI , also the opposite velocity v* = —v occurs with the 
same probability (under the tacit assumption of an analo- 
gous equal a priori likelihood of the initial condition x% (0) 
and x*(0)). Since, as mentioned above, different initial 
conditions indeed may lead to different velocities, it fol- 
lows that our symmetric system may give rise to solutions 
which break the symmetry of the system by exhibiting 
a finite velocity in cither one or the other direction. As 
usual (see above), we will refer to such a case as sponta- 
neous symmetry breaking and to the concomitant dimcr 
motion as spontaneous transport. 

Turning to T > (and arbitrary F), the above men- 
tioned fact that the velocity v from (fT3"|) is independent 
of the initial conditions and the realization of the noise 
results in the basic symmetry property that F M> — F 
implies v 1— > — v. 

In particular, in the unbiased case F — we recover 
v = for any T > 0. Regarding finite bias forces F, one 
could naively expect that the velocity (TIB"]) now should 
take on a finite value and that its sign should always 
correspond to that of F. The quite unexpected contrary 
behavior, namely a velocity v opposite to the applied 
force, is termed absolute negative mobility [54|. 

The main aim of our paper is to demonstrate that the 
dynamics (TT|) indeed may give rise to spontaneous trans- 
port and to absolute negative mobility, and that there 
exists in fact a close connection between these two quite 
astonishing effects. A comparison with related previous 
works is provided in Sec. X. 

IV. NO-GO THEOREMS FOR CONVEX 
INTERACTION POTENTIALS 

In this section we demonstrate that spontaneous trans- 
port and absolute negative mobility (as specified in the 
previous section) are ruled out for any dynamics p])-(|4| 
with a convex interaction potential W(y), i.e. 

W"(y) > for all y. (16) 

In order to exclude spontaneous transport, it is suf- 
ficient to demonstrate that the velocity (|13p in the de- 
terministic limit (T = 0) is unique (independent of the 
initial conditions). Proceeding by way of a proof by con- 
tradiction, we assume that there exist two solutions of 
dU)-OJ> , called Xi(t) and x'^t), with different velocities v 
and v' according to (TT3")) . say v 1 > v without loss of gener- 
ality. Since for any integer n, also Xi(t)+nL solves (jT|)-(j4j) 
and yields the same velocity v in (|13[) . we may without 
loss of generality assume that the initial conditions sat- 
isfy xi(0), £2(0) > x[(0), x' 2 (0). Since v' > v, the parti- 
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cles xj(t) must "catch up" with the particles Xi(t) in the 
course of time. Without loss of generality, the indices i 
can be chosen so that it is x[(t) which catches up with 
X\{t) for the first time at some t = t* , i.e. 

x[(n = xi(f) (17) 

x' 2 (t*) < x 2 (t*) (18) 
±i(**) > x x (t*) . (19) 

The equality sign in (|18[) can be excluded, since otherwise 
(|T7)) and the dynamics ©-flU) (with T = 0) would imply 
that a^(t) = for all i, and thus i>' = v. With j/ := 
x\{t*) - x 2 {t*) and y' := x'^t*) - x' 2 {t*) it follows that 

y'>y, (20) 

and by introducing jTT} into (H}-® (with T = 0) we find 

±i(**)-±i(**) - + ■ (2i) 

In view of (|2"0)) and (|16p we thus can infer that x'^t*) — 
xi(t*) < 0, in contradiction to (fT9")l . This concludes our 
proof that spontaneous transport is excluded. 

Next wc analogously proof by contradiction that (fl6|) 
rules out absolute negative mobility for the dynamics ([1])- 
(|4]) with T > 0. We thus consider a solution Xi(f) of 
(HI)— d3J) with velocity v according to (|T3"|) and a second so- 
lution x\ (t) with F' instead of F in (JTJ) and a resulting ve- 
locity i/. Assuming absolute negative mobility, it follows 
that two such solutions must exists with the property 
that F' < F and v' > v. In particular, we can assume 
the same realization of the noise & (£) in (TT]) for X{ (t) and 
for x'^t), since - as mentioned in Sec. Ill - the velocity 
([13]) is independent of the realization of the noise. As be- 
fore, Xi(0), x 2 (0) > x[(0), x' 2 (0) can be taken for granted 
without loss of generality and v' > v then implies the ex- 
istence of a time point t* with the properties (|T7|) - (rT9l . 
With y := xx(t*) - x 2 (t*) and y' := x[(t*) - x' 2 (t*) it 
follows that 

y'>y, (22) 

but, unlike before, y' = y cannot be excluded now. From 
(HJ-Q we can infer that 

x'^t*) - Xl (t*) = -W'(y') + W'{y)+F' -F .(23) 

Combining gjo§ with (J22J> yields -W'(y') + W'(y) < 
and with F' < F it follows that x[(t*) - ±i(t*) < 0, in 
contradiction to (TT5]). We thus have proven that absolute 
negative mobility is impossible. 

We remark that the above arguments can be readily 
generalized to more than two and/or non- identical parti- 
cles with convex nearest-neighbor interaction potentials. 
In particular, in the context of the Frenkel-Kontorova 
model, the so-called Middelton's no passing rule has been 
proven along similar lines, see [55l459| and further refer- 
ences therein. 



V. NON-CONVEX INTERACTION 
POTENTIALS 

Non-convex interaction potentials occur naturally in 
various situations: First, this is the case for many com- 
mon models of basic molecular interactions [J-Q , such as 
Lennard- Jones [67[ (cf. ©) or Morse potentials fsl. |34|. 
Second, less basic but still bistable molecular interaction 
potentials arise, e.g., in the modeling of random walkers 
in motor proteins |9h181 l6r| , and have also been consid- 
ered in dimer ratchet models [H, H3, 0, S3 ■ Third, more 
complex dimers may have several stable configurations, 
or the dissociation process may involve additional (fast) 
internal degrees of freedom that lead to a non-convex ef- 
fective interaction potential, as considered theoretically 
e.g. in [53}. Experimentally, e.g. in [8(j a water droplet 
was positioned on a periodic surface and its shape (in- 
ternal degree of freedom) was periodically driven by ex- 
ternally applied electric fields. Lastly, a spatially one 
dimensional model may serve as an approximation for 
the full three-dimensional dynamics of the dimer. The 
main effects of the other degrees of freedom are incor- 
porated into the "effective" interaction forces between 
the monomers, resulting in a more complex, non-convex 
structure of the potential (cf. Sec. IX). 

VI. SPONTANEOUS SYMMETRY BREAKING 

We consider the dynamics (HJ)— (I3J) without bias (F = 0) 
in the deterministic limit (T — 0). We furthermore focus 
on non-convex interaction potentials W(y), for instance 
(lU), ([9]). Is it possible that a solution Xi(t) spontaneously 
breaks the spatial symmetry of the system and exhibits a 
finite velocity v according to Eq. (fT"3"|) ? While there seems 
to be no a priory reasons why not, the actual existence of 
such solutions appears to us quite astonishing neverthe- 
less. In accordance with this intuitive expectation, their 
existence turns out to be indeed restricted to relatively 
small parameter regions in the deep non-linear regime, 
thus making any further progress by analytical means 
quite hopeless. Therefore, we focus on the discussion of 
our extensive numerical results in the first place and only 
a posteriori provide a simplified intuitive picture of what 
is going on. 

A. Numerical example 

Figure 3 presents typical results for the quartic poten- 
tial from ([9]) (see also Fig. 2b). For each depicted point 
("pixel") in Fig. 3, we have calculated the velocity ac- 
cording to (fT3"|) for 20 trajectories of length 500 r with 
randomly chosen initial conditions. To eliminate "tran- 
sient" effects, we furthermore excluded an initial time- 
interval of length 150 r for each trajectory contained in 
the average in (fT3")) . The black "islands" indicate pa- 
rameter regions where all initial conditions which we nu- 
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merically tried out ended up on periodic attractors ex- 
hibiting spontaneous symmetry breaking, i.e. a finite net 
velocity of the rational form v = ±n/m with positive in- 
tegers ri, m in the specific dimensionless units from (|1 3|) . 
We note that spontaneous symmetry breaking may occur 
also with other types of attractors, e.g. chaotic attractors 
which e merg e from the periodic ones by period doubling 
cascades [63[ and thus still carry a rational net velocity of 
the form v = ±n/m. However, they are not visible in Fig. 
3, as their parameter regions are extremely small. Fur- 
thermore, we found that the black islands in Fig. 3 are 
largely dominated by the velocities v = ±1. The light 
grey "sea" around the black islands indicates parame- 
ter regions without spontaneous symmetry breaking, i.e. 
all probed initial conditions exhibited a net velocity of 
V = (in fact, every trajectory x(t) remained bounded, 
and, in particular, did not exhibit deterministic diffu- 
sion). The quite small remaining regions in mid-dark 
grey indicate coexistence of solutions with v = ±n/m 
(like in the black islands) and with v = (like in the 
light grey sea). The tiny (almost invisible) white regions 
and single pixels in Fig. 3 indicate that the dimer tra- 
jectory was neither bounded nor did it reach a periodic 
attractor within the maximal numerical simulation time 
of 1000 r. Since those parameter regions are so small, we 
will not consider them any further. 

The negative parameter regimes (omitted in Fig. 3) are 
cither unphysical (e.g. c in (O), follow by symmetry (e.g. 
A and u) in pT]) ). or by periodicity (e.g. A in ©). For 
all the remaining (non-negative) parameters outside the 
range plotted in Fig. 3, either only solutions with v = 
were obtained, or - in the case of large A- values - the con- 
tinuation is qualitatively quite obvious. By changing the 
values of those parameters which are held fixed in each 
of the three plots in Fig. 3, also the shape and positions 
of the "black islands" change considerably. E.g., we have 
found that for every given 0.001 < uj < 2 it is possible 
to adapt the remaining parameters so that spontaneous 
symmetry breaking is realized (even smaller ui- values be- 
come numerically very expensive and have therefore not 
been investigated any further). From this viewpoint, the 
range of symmetry breaking parameters is actually not 
so small as a first glance on Fig. 3 may suggest. 



B. Basic mechanism 

In order to illustrate the characteristic features of the 
dimer motion when the symmetry is spontaneously bro- 
ken, we focus, as a typical example, on the parameter 
values of the driving indicated by the "+" symbol in 
Fig. 3a and on the specific interaction potential W(y) 
depicted in Fig. 2b. For these parameter values, two 
periodic attractors with net velocity v = ±1 coexist, 
which break the symmetry of the underlying equations 
of motion (|lll) , (fT2"T) . Figure 4 collects 9 snapshots of the 
associated dynamical evolution of the dimer during one 
driving period r on the attractor with v = +1. At the 
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FIG. 3: (Color online) Velocity (|13p obtained by numerically 
solving |TT]), (H2j) with F = 0, T = 0, and © for various 
parameter values, (a) Variable driving amplitude A and fre- 
quency u> with fixed parametes c = 0.428, A = 0.26L (these 
"odd" values are chosen on purpose to avoid all "acciden- 
tal symmetries" or other sorts of "special cases"). The two 
crosses "+" and "x" indicate the specific parameter values 
adopted later in Fig. 4 and Fig. 6, respectively, (b) Variable 
driving amplitude A and interaction potential parameter c 
with fixed parameters ui — 0.064, A = 0.26 L. (c) Variable 
driving amplitude A and interaction potential parameter A 
with fixed parameters ui — 0.064, c = 0.428. Black: Finite 
velocities v = ±n/m with positive integers n and m. Light 
grey: Velocities v = 0. Intermediate grey: Coexistence of fi- 
nite velocities v = ±n/m and v — 0. For further details see 
main text. 



beginning of the cycle, the dimer extension is close to its 
equilibrium value of y = L/2 = tt (see also Fig. 2b and 
Eq. (pjj)- In this configuration, the forces due to the ex- 
ternal periodic potential acting on the two particles are 
almost identical in modulus but have opposite signs, so 
that the center of mass motion of the dimer is practically 
not affected by the periodic potential (cf. Eq. (fTTj) with 
y rj L/2 = 7r), whereas its impact on the relative coordi- 
nate y is maximal (cf. Eq. (fT2|l). As a consequence, the 
dimer as a whole essentially follows the time-dependent 
driving force, while it is slightly stretched and compressed 
by the periodic potential landscape. This type of be- 
havior can be observed in Fig. 4 during the first half 
of the driving period r and results in a displacement of 
the dimer's center by a little more than one spatial pe- 
riod in the positive x direction. The deviations from the 
equilibrium length y — L/2 are stabilized very quickly 
due to the convex character of the interaction potential 
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FIG. 4: (Color online) Time evolution of a typical spon- 
taneous symmetry breaking periodic attractor of (jlll) . (|12[) 
with interaction potential @ and the specific parameter val- 
ues indicated by the "+" symbol in Fig. 3a, i.e. uj = 0.175, 
A — 0.885 and all remaining parameters as in Fig. 3a. The 
black dots indicate xi(t) (right dot) and X2(t) (left dot), the 
arrows x\(t) and X2(t), and the black lines the tilted periodic 
potential U(x) — (x — L/2)f(t). The dark grey "insets" in- 
dicate the interaction potential ([9| (grey line) along with the 
instantaneous dimer length y(t) (grey dot). 



W(y) for y w L/2. However, the non-convexity of W(y) 
for j/-values around y = A = 0.26L allows the dimer to 
be eventually compressed to about L/4 in length, and to 
maintain this "metastable" state for almost the complete 
second half of the driving period. During this phase, the 
dimer is trapped motionless in a minimum of the rock- 
ing potential U{x) — xf(t). At about 6r/8, when the 
left potential barrier of the minimum is nearly balanced 
by the external driving (although the minimum actu- 
ally never disappears because the driving amplitude A 
is smaller than unity), the dimer starts to relaxate back 
to its equilibrium length y = L/2 by way of "pushing 
the left monomer over the barrier" . In this manner, the 
dimer reaches its initial configuration, but with the cen- 
ter being shifted by n — 1 spatial period to the right after 
m = 1 cycle of the driving, resulting in v = +1. 

From these observations we get an intuitive idea, in 
addition to the formal arguments of Sec. IV, of why non- 
convexities in the interaction potential arc necessary for 
spontaneous symmetry breaking. They promote the ex- 
istence of a "metastable" dimer state with dynamical be- 
havior different from that of the equilibrium configura- 
tion. The dimer motion in the periodic potential land- 
scape under the action of the external driving induces 
periodic transitions between that "metastable" state and 
the equilibrium configuration, associated with distinct 



phases in the dimer motion and thus resulting in an over- 
all symmetry breaking behavior. The above example has 
been chosen to most clearly illustrate this mechanism, 
as the specific parameter values give rise to an equi- 
librium configuration with practically "free" motion of 
the dimer's center, whereas in its "metastable" state, the 
dimer is trapped motionless by the potential landscape. 
Such a clear-cut distinction between the two states is, 
however, not necessarily required. Figure 3 shows that 
spontaneous symmetry breaking occurs in quite large pa- 
rameter regions, where the "metastable" state is "suffi- 
ciently different" from the equilibrium configuration of 
the dimer. 

From our discussion so far, it is evident that average 
dimer velocities larger in modulus than \v\ = 1 should 
be observable. For instance, for slower driving (smaller 
w) the dimer may be able to travel several spatial pe- 
riods during the first half of the driving cycle, while it 
gets trapped during the second half of the cycle. In- 
deed, dimer velocities with v = ±2 are present in the 
black regions of Fig. 3a for small ui and small A val- 
ues (in Fig. 3, the information on the actual value of 
non-vanishing average velocities in the black regions is 
not included though). Likewise, rational velocities with 
\v\ < 1 can, e.g., be found for fast driving at the edges of 
the black regions towards large ui and A values. 

While these considerations are helpful to gain an intu- 
itive picture of what is going on, they obviously do not 
prove that a specific symmetry breaking solution is dy- 
namically stable (an attractor) and that, besides its twin 
brother with opposite velocity, no other stable solutions 
exist. Any such "more rigorous" reasoning would have 
to depend on the precise quantitative values of all model 
parameters and thus seems almost impossible. 



C. Generalizations 

We finally point out that the specific shape of the 
interaction potential is not crucial for the existence of 
spontaneous symmetry breaking, as long as it is non- 
convex. For instance, we have found spontaneous sym- 
metry breaking attractors for various modifications of 
©, including the obvious situation where it has a sec- 
ond minimum, and also the "extreme" case of an "al- 
most convex" potential with a curvature that is positive 
everywhere except near the equilibrium length y ss L/2. 
Moreover, spontaneous symmetry breaking is also sup- 
ported by completely different non-convex interaction po- 
tentials. 

As a specific example, we now turn to our numeri- 
cal results for the Lennard- Jones potential (J8J), see also 
Fig. 2a. For large distances y = x\ — X2 of the monomers, 
their interaction becomes negligibly small. As expected, 
we found numerically that in this case each monomer 
exhibits a periodic back-and-forth motion in the long 
time limit and thus the velocity v from Q13p vanishes. 
Since these solutions are of little interest, we hence- 
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FIG. 5: Same as Fig. 3a but for the Lennard- Jones interaction 
potential ([8]) with A = 1.25 L. Only attractors that can be 
reached by dimers initially close to their equilibrium lengths 
are shown, see text. 




forth use only initial conditions such that the dimer 
length is initially close to the equilibrium length A of 
the potential ©. More precisely, we only considered 
y(0) € [0.9 • A, 1.3 ■ A] and observed that the distance 
y{t) always remained reasonably small to consider the 
dimer as "not dissociated" (a rigorous dissociation cri- 
terion for the potential © seems hard to define). The 
so obtained numerical results are qualitatively very sim- 
ilar to those for the quartic potential ([9]), although the 
quantitative details arc different. As an example, Fig. 5 
depicts our findings for the Lennard- Jones potential ([8]) 
corresponding to those from Fig. 3a for the quartic po- 
tential ©. In particular, we recover spontaneous sym- 
metry breaking within a quite notable parameter range. 
The first main differences compared to the quartic poten- 
tial is that we did not find coexistence of solutions with 
v = and u ^ 0. The second main difference is that all 
the spontaneous symmetry breaking solutions that we 
found exhibited the specific velocities v = ±1 for all the 
parameters that we tested. 



VII. GIANT DIFFUSION 

We still focus on the unbiased (F — 0) dynamics ([]}- 
(0, but now consider finite noise strengths (T > 0), 
implying v = for all model parameters and all initial 
conditions (see Sec. III). The reason is that any attrac- 
tor of the deterministic dynamics - in particular those 
giving rise to spontaneous symmetry breaking - now be- 
comes mctastable due to noise induced transitions be- 
tween them. While the velocity (fT3")) is thus trivial, the 
quantity of foremost interest is now the diffusion constant 

®. 

For large noise (T > 1) in ©-©, the effects of 
the "substrate potential" U (x) become negligible, hence 
the diffusion asymptotically approaches the value for 
the "free" thermal diffusion of the dimer, namely D = 
k B T/2r]. 

Much more interesting is the opposite asymptotics of 
small T, specifically in the case of spontaneously bro- 
ken symmetry without coexisting deterministic solutions 
with v = 0, as exemplified by the black islands in Figs. 



FIG. 6: (Color online) Arrhenius plot (logarithm of D versus 
1/T) of the diffusion coefficient (|14[) by numerically solving 
the unbiased (F = 0) dynamics (|lip . (|12|l with quartic inter- 
action potential ([9]) and the specific parameter values indi- 
cated by the "x" symbol in Fig. 3a, i.e. u = 0.11, a — 0.95 
and all remaining parameters as in Fig. 3a. 



3 and 5. Focusing on this case, a solution Xi(t) of the 
noisy dynamics <HJ)-((1]) closely follows most of the time 
cither a deterministic solution exhibiting v > or its twin 
brother with v < 0, switching between them every now 
and then by way of noise activated transitions. For suf- 
ficiently weak noise, these transitions become very rare 
events and are captured very well by a rate process with 
a certain transition rate k. A similar calculation as in 
Ref. [8l| then yields for the diffusion coefficient (ll~4l the 
result 



D = 



2r 2 k 



(24) 



The transition rate itself is of the usual Arrhenius-form 



k = v a cxp{-A0/T} 



(25) 



with a T-independent generalized activation energy Acj) 
[73l ] and an "attempt frequency" Vq with a rather weak 
T-dependence (at most algebraic). Combining (I24p and 
(1231) . we thus arrive at the analytical prediction [82| of 
an extremely strongly diverging diffusion coefficient for 
asymptotically weak noise strengths T (essential singu- 
larity). The numerical findings in Fig. 6 confirm this 
prediction very well. 



VIII. ABSOLUTE NEGATIVE MOBILITY 

As before, we consider the dynamics (JTJ)— with T > 
0, but now also in the presence of a non- vanishing static 
bias F. For symmetry reasons (see Sec. Ill) we hence- 
forth focus on F > 0. Our main question is: how does 
the unperturbed velocity (v = for F = 0) respond to a 
perturbation F > 0, and, in particular, is it possible to 
obtain absolute negative mobility, i.e. v < ? 



To be specific, we revisit the deterministic (T — 0) and 
unbiased (F = 0) dynamics (fTTj) , (fT2"|) for the quartic po- 
tential (|9|), resulting in the velocity diagram from Fig. 3a. 
We recall that every black island in this figure indicates 
the coexistence of two symmetric deterministic attrac- 
tors, one with v > and the other with v < 0, and that 
all solutions (for arbitrary initial conditions) approach 
one of these two attractors in the long time limit, i.e. all 
solutions exhibit either v > or v < 0. As soon as F > 0, 
the symmetry of the system is broken and hence the de- 
generacy of the u / solutions lifted. Accordingly, by 
gradually increasing F, the "island" in parameter space 
describing solutions with v > will start to differ more 
and more from the v < "island" . Every subset of the 
latter island with v < which is located "outside" of 
the v > island is clearly a quite promising candidate 
for absolute negative mobility at finite (sufficiently small) 
T. Fig. 7a confirms that such subsets indeed exist and 
Figs. 7b, c show that they indeed give rise to absolute 
negative mobility. 

By comparing Fig. 7a with Fig. 3a we see that the de- 
generacy of the v > and v < solutions is lifted by dis- 
placements of the red v > "islands" toward larger A and 
oj values and of the cyan v < "islands" toward smaller 
A and oj values. While the impact of a change in the driv- 
ing amplitude A on a specific solution and its stability is 
hard to grasp intuitively [83[ , the existence of a "purely 
cyan stripe" outside the border of any "red island" to- 
ward smaller frequencies oj can be understood as follows: 
In the unbiased case (F = 0) from Fig. 3a, the lower bor- 
der of an island indicates that solutions with v ^ will 
cease to exist (become unstable) upon further decreasing 
the driving frequency oj. Roughly speaking, the driving 
thus becomes "too slow" compared to the "eigenspeed" of 
the solutions, i.e. they are unable to synchronize with the 
imposed temporal periodicity of the driving. Moreover, 
if a small bias F > is switched on, it is quite plausible 
to expect that the intrinsic "eigenspeed" will decrease for 
solutions traveling "uphill" (against the force F), and in- 
crease for those traveling downhill. As a consequence, at 
the lower border of an island in Fig. 3a, those solutions 
with v < should survive and those with v > disap- 
pear when a small F > is switched on, in agreement 
with Fig. 7a. 

Upon further increasing F, one expects that the cyan 
and red islands in Fig. 7a change not only their posi- 
tions, but also their shapes and in particular their areas, 
and that for sufficiently large F, the cyan islands will 
completely disappear. Likewise, with increasing noise 
strength T, one expects that the original deterministic so- 
lutions become less and less relevant. Indeed, for T — > oo 
the effects of the substrate potential U (x) in (|4]) become 
negligible, yielding for the velocity (|T3| the trivial result 
v = F/rj. These quite plausible qualitative features are 
indeed confirmed by our numerical results for a represen- 
tative example in Fig. 7c. 

Turning to the Lennard- Jones potential (|8|). it is quite 
plausible and can be confirmed by more detailed calcu- 




F 



FIG. 7: (Color) (a): Same as Fig. 3a, but for a finite bias 
F = 0.01. Red: Velocities v > 0. Cyan: Velocities v < 0. Red 
and cyan overlapping: Coexistence of solutions with opposite 
average velocities. Grey: Velocities v = 0. Coexistence of 
v and v = is resticted to negligibly small parameter 
regions and is not shown in this figure, (b): Same as in (a) 
but for a finite noise strength T — 0.005. The value of the 
velocity v from (|13jl is independent of the initial conditions 
and indicated by the color code. In particular, blue regions 
correspond to absolute negative mobility (v < for F > 0). 
(c) Dependence of the velocity v on the bias F and the noise 
strength T for fixed driving parameters as indiciated by the 
black "+" in (a) and (b). 



lations that any solution of ([TJ with T > results in a 
dissociation of the dimer in the long run, and that this 
implies v > for all F > 0. However, for small T and F, 
the lifetime of the dimer is still very large, and restricting 
the time-average of the velocity in (fTU| accordingly, the 
same qualitative features as in Fig. 7 are recovered (not 
shown). Fig. 8 exemplifies a typical solution exhibiting 
absolute negative mobility during the extremely long ini- 
tial time-period before the dimer dissociates. The very 
long life-time of the dimer can be understood intuitively 
as follows: As a relict of the underlying (dynamically sta- 
ble) symmetry-breaking attractor, the dimer length never 
increases notably beyond 2 during the motion shown in 
Fig. 8. The binding energy of the dimer at that length 
is still ©(lCP 1 ) (see Fig. 2b) and thus much larger than 
the thermal energy corresponding to T = 0.001, so that 
noise-induced dissociation events are extremely rare. 
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t/T 

FIG. 8: (Color online) (a) Center of mass x(t) from (fT0|) 
for the dynamics JTTJ), (|12|) with Lennard- Jones interaction 
potential ([9]). Remaining parameters: A = 1, u = 0.06, 
T = 0.001, F = 0.01. (b) Dimer length y(t) at the time 
points t = n/10 for n = 0, 1, 20000. 



IX. OUTLOOK ON HIGHER DIMENSIONAL 
MODELS 

In this section we address the motion of two coupled 
particles in 2- or 3-dimensional periodic potentials, e.g. 
on a crystal surface, in a lattice of optical or magnetic 
traps etc., see also Sec. I. A particularly interesting spe- 
cial case arises for systems with a pronounced anisotropy, 
causing a strong preferential direction ("easy axis") for 
the dimer motion and orientation, as exemplified with 
Fig. 1. 

Focusing on this case, the one-dimensional model ([1])- 
(|4|) still describes such systems reasonably well, espe- 
cially with the following modification: The main effect of 
the omitted additional dimensions is that the dimer may 
"turn around" by 180°, albeit with a very low probabil- 
ity per time unit (due to the strong anisotropy). This 
feature can be effectively incorporated into the model 
([T])-((4]) by choosing an interaction potential W(y) which 
exhibits two minima, one at a positive and one at a neg- 
ative y- value, separated by a "high" energy barrier at 
y = (to guarantee "rare" transitions). In other words, 
a sign change of x\ — X2 = y (cf. Eq. (flO|) ) does not 
mean that the two particle pass "through" each other, 
but rather that they perform a rotation. In contrast to 
our convention at the end of Sec. II, we furthermore re- 
strict ourselves to "physically realistic" interactions, de- 
pending only on the absolute particle distance and thus 
satisfying 

W(-y) = W(y) . (26) 

In other words, W{y) exhibits a symmetric double well 
structure and thus is always a non-convex potential. 

Numerically, we found that effectively one-dimensional 
models of the above type reproduce all the main ef- 
fects from Sees. VI- VIII. Essentially, the reason is that 
once the solutions of a model as specified in Sec. II 
are known, the solutions of the same model but with 



W(y) i — ^ W(y + L) follow immediately and the resulting 
velocities and diffusion coefficients remain unchanged. 
Hence, our present double well potential W(y) with a 
barrier at y = can be readily shifted by any integer 
multiple of L "to the right" until it becomes similar to 
the one of the previously considered interaction poten- 
tials (see Sec. VI.C). 

More generally, also dimers consisting of non-identical 
particles can be approximated by an effectively one- 
dimensional model analogous to ([I])- (HI, except that the 
friction coefficients r; in (JT)) may now be different for the 
two particles, and similarly for the "coupling strengths" 
to the forces F, f(t) and the substrate potentials U(x). 
As a consequence, one can show that the system satis- 
fies only a weaker symmetry property than the one from 
Sec. III. Accordingly, a spontaneous breaking of this sym- 
metry is in many cases rather trivial: Namely, the dimer 
exhibits a standard ratchet effect since the two particles 
arc different and never actually exchange their positions 
during their time-evolution [3J, H(| [47}. Considerably 
less trivial at first glance is the occurrence of giant dif- 
fusion and absolute negative mobility in such systems. 
But given our previous discussions in Sees. VII and VIII, 
these effects are, of course, immediate consequences of 
the spontaneous symmetry breaking. 



X. CONCLUSIONS 

We considered a very simple "minimal" model arising 
in a large variety of different situations (see Sec. I): Two 
coupled ovcrdamped Browian particles moving along a 
periodic substrate potential in one or more dimensions. 
In spite of the fact that the particles are identical and 
the substrate spatially symmetric, a sinusoidal external 
driving may lead to spontaneous symmetry breaking in 
the form of a permanent directed motion of the dimer 
(spontaneous transport, see Sec. V) in the deterministic 
limit (T = 0). Thermal noise (T > 0) restores crgod- 
icity and thus zero net velocity, but entails arbitrarily 
fast diffusion of the dimer for asymptotically small tem- 
peratures T (Sec. VII) Moreover, for small but finite T 
the dimer responds to a static bias force with motion op- 
posite to that force (absolute negative mobility, see Sec. 
VIII). As proven in Sec. IV, the key requirement for all 
these effects is a non-convex interaction potential of the 
two particles (Sec. V). 

Directed transport induced by spontaneous symmetry 
breaking as a collective effect of infinitely many particles 
far from equilibrium has been observed before, e.g., in 
Refs. [l||, [6J, [6^, H|[ . These systems and also the basic 
physical mechanism are quite different from our present 
case (Sec. VLB). Spontaneous symmetry breaking in- 
duced transport for periodically driven, underdamped 
(finite inertia) systems [69|-[72j is, mathematically speak- 
ing, closer related to our case (both types of dynamics 
are two-dimensional, dissipative, and non-autonomous), 
while the underlying physical systems (Sec. I) are en- 
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tircly different. 

Diffusion enhancement phenomena are well established 
already for a single overdamped particle under the influ- 
ence of an unbiased, temporally periodic driving force 
[60l - [62l ] . But a diverging diffusion coefficient for asymp- 
totically weak noise has been reported for the first time 
only very recently in 82j. In particular, we remark that, 
unlike in (60T - [62| . in our present case, the parameter space 
of diffusion enhancement does not approach a set of mea- 
sure zero for asymptotically weak noise. 

Absolute negative mobility has been reported in a con- 
siderable variety of experimental systems and theoretical 
models by now, see j54j for a recent review. Within the 
realm of classical models (no quantum effects), it has 
been observed for the first time as a collective effect of in- 
finitely many interacting particles under far from equilib- 
rium conditions (74[. Later, this model was successfully 
reduced down to 3 interacting particles |75j . Although 
this model is very different from ours, the present results 
for two interacting particles may be seen as the next and 
ultimate reduction step, since, as demonstrated in [7l| , a 
single overdamped particle in one dimension cannot ex- 
hibit absolute negative mobility. 

A somewhat similar stochastic dynamics of two cou- 
pled Josephson junctions exhibiting absolute negative 
mobility has recently been studied in Ref. [29|. This 
model is mathematically equivalent to the overdamped 
dynamics of two coupled particles. However, it exhibits 
certain features which are very different from our present 
model: (i) The two particles are not identical in all re- 
spects; (ii) The interaction forces cannot be derived from 
an interaction potential. Moreover, they are spatially 
periodic, i.e. not approaching zero for large particle dis- 
tances. Besides those differences in the considered model, 



also the basic mechanism generating absolute negative 
mobility is very different. 

Further, mathematically somewhat similar but phys- 
ically completely different systems exhibiting absolute 
negative mobility are (i) a single underdamped (finite 
inertia) particle in one dimension [tTL I72I l7tSj . (ii) an 
overdamped particle on a one-dimensional "meandering" 
state space [Hil. l77j (iii) an overdamped particle in 2 di- 
mensions [78j , and (iv) a single particle in one dimension 
with some kind of "internal degree of freedom" [79| . The 
recent work [52j may be considered as an extension of 
the above mentioned system class (i) to the case of two 
coupled underdamped particles. 

In the context of nano-friction experiments b y a tomic 
force microscope (AFM), it has been reported 84J that 
an external time-periodic forcing may result in unmea- 
surably small sliding friction forces (so-called superlu- 
bricity). Given that dimer models are well established 
to describe such systems (22T-|25j , our present results pre- 
dict the possibility of an even more astonishing effect, 
namely sliding friction with negative sign, i.e. opposite 
to the direction of the pulling velocity. 

For colloidal particles in magnetic or optical traps [2(| 
l2~l| , an effective "bias" force is often generated by moving 
the traps with respect to the ambient fluid and then going 
over into the comoving reference frame. Along the lines 
of our present work, it might be possible to make colloidal 
doublets [l9| move oppositely to the moving (and usually 
"entraining" ) traps. 
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